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ABSTRACT Investigations on diffusion in systems with memory [1] have established a 
hierarchical connection between mixing, ergodicity, and the fluctuation-dissipation theorem 
(FDT). This hierarchy means that ergodicity is a necessary condition for the validity of the 
FDT, and mixing is a necessary condition for ergodicity. In this work, we compare those 
results with recent investigations using the Lee recurrence relations method [2-4]. Lee 
shows that ergodicity is violated in the dynamics of the electron gas [4]. This reinforces 
both works and implies that the results of [1] are more general than the framework in 
which they were obtained. Some applications to slow relaxation phenomena are discussed. 

Introduction - More than a hundred years after its formulation by Boltzmann, the ergodic 
hypothesis (EH) still drives the attention of the mathematics and physics community. Many 
situations have been found in which the EH is broken [1, 3-5]. On the other hand, the 
mixing condition (MC), despite being omnipresent in relaxation processes, has seldomly 
been the object of investigation [1,6]. Recent studies on anomalous diffusion [1,7] have 
established a strong hierarchy, which in growing generality is: the FDT, EH, and the MC. 
Since the FDT is a less general "theorem", we shall focus only on the EH and the MC. 

Kubo [8] realized that doing a time average on correlation functions would be more realistic 
than performing those averages on the variables themselves. He then established an ergodic 
condition in the framework of linear-response theory. That is, if the zero-frequency limit of 
a dynamical susceptibility is equal to its static counterpart, the system is ergodic. However, 
the problem of solving the dynamical equations of motion for many-body systems remains 
a daunting one at best, consequently the validation of the EH via Kubo's method has been 
done only for a few systems. 

With the development of the recurrence relations method by Lee [2], it became possible to 
obtain general solutions for the Heisenberg equations of motion. In this way it is possible 
to verify the validity of the EH [3, 4] .We have avoided this problem using the asymptotic 
limit of the correlation function, which allows us to verify the EH without a full solution of 
the dynamical equations of motion. For that, we just need the memory, which is our starting 
point. In this work we try to narrow the difference between the methods, i.e. we show that 
both methods yield the same result. 

Ergodicity - Recently, Lee [4] proposed that the validation of the EH is subject to the 
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condition 



< r = / R(t)dt < oo, (1) 
J o 

where r is a relaxation time we introduce here for further use, R(t) is the renormalized 
correlation function for the dynamical stochastic variable A(t). It is given by 

_ (A(t)A(0)} 
m ~ (AW) ( 

and here (A(t)A(0i)) is the ensemble average. The finite value of the integral is the Lee 
condition [4] for the EH be valid. Note that a necessary condition for the validity of Eq. Q 
is that the mixing condition (MC) holds, i.e. 

lim R(t) = 0, (3) 

t— >oc 

otherwise the integral would diverge. Consequently the MC is necessary for the EH, see 
Costa et al. [1]. The condition imposed by Lee is sufficient for the validity of the EH. 
However, there are situations where it may be too restrictive and this condition can be 
weakened. We shall prove here that in anomalous diffusion the MC is a sufficient condition. 

Anomalous diffusion can be well described by a generalized Langevin equation (GLE) of 
the Mori form [1,7]. From the GLE it is possible to obtain 



dR(t) 



[ U(t-t')R(t')dt', (4) 
J o 



dt 

where IL(t) is the memory. The Laplace transform of this equation yields 

R(?) = ^TT> (5) 

z + U(z) 

where the tilde indicates Laplace transforms. Consider 

lim R(t) = lim zR(z) = lim 1 , (6) 



i^oo z^O z^Q 



U(z) 



where we have used the final value theorem [9]. Note that the inverse Laplace transform 
can be obtained analytically only in a few cases. However, since we know ti{z) the limit 
can be obtained even without an explicit solution for R(t). For 

U(z - 0) ~ z v , (7) 



the MC is not valid for v > 1. For diffusive processes with the mean-square displacement 
given by (x 2 (t)} oc t a , the exponent is [1, 7] a = v + 1. For v = — 1 we have a = 0, 
i.e. no diffusive process. This case needs special analysis: let v = —1, i.e. T(z) = K/z. 
The inverse Laplace transform gives F(t) = K. A constant value will produce a harmonic 
oscillator term which localizes the particle, i.e. no diffusion, in agreement with a = v + 1. 
From Eq. (0, we get that R(t) = cos(VKt). We can see that this correlation function 
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has an infinite relaxation time, and violates the MC and EH. Consequently, mixing and 
ergodicity in diffusive processes will hold only for —1 < v < 1, which correspond to 
< a < 2. For v = (normal diffusion), we can define a finite relaxation time by 
Eq. 0. Strictly speaking, one can define a relaxation time only for normal diffusion. In 
that case, with broad-band noise, the correlation function relaxes as exp(— t/r) [10]. We 
shall analyze here the extreme cases where the MC condition is violated. 

Slow Relaxation Phenomena - We shall restrict ourselves to diffusive phenomena in order 
to compare our results with those of Lee. The majority of systems that violate ergodicity 
present some form of slow dynamics. For example Mukamel et al. [5] show that Ising sys- 
tems with long-range interactions exhibit a relaxation time which diverges logarithmically 
with system size. 

Another example of violation of the MC, which again implies the violation of the EH, is 
the ballistic motion, which we shall discuss in the next section. Let us now address the 
superdiffusive motion i.e. < v < 1. 

As we have mentioned before, relaxation times exist only for normal diffusion, However, 
for superdiffusion, Eq. © in the limit given by Eq. (0, yields a Mittag-Leffler function 
of the form E\- u {— (t/T v ) l ~ u ), which displays a transient behavior from a stretched ex- 
ponential exp — (t/r u ) 1 ^ u to a power law {t/T v -\) u ~ l . In order to have an idea of the 

magnitude of r, we need to know the behaviour of n(z) for small z. This depends on the 
nature of the noise in the stochastic process [6] For a noise density of the form 
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, if LO < UJ£> 

(8) 

0, otherwise, 



with ujd as a Debye cutoff frequency, it is possible to compute the memory [1,7] using the 
Laplace transform, and its low z limit, to obtain the coefficient of Eq. Q. The transient 
time is t u = tq [(^/uidY sec ((vit)/2)] 1 ^ u ~ 1 \ In Fig.[Qwe plot the transient time t v as a 
function of v for several values of ujd/i- For normal diffusion, to = l/7> it i s equivalent 
to the relaxation time. Notice that the maximum increases with ujdH- F° r broadband noise 
ud/i S> 1, the transient time becomes very large as v approaches 1. 

Ballistic Motion - In nature, normal diffusion and subdiffusion are prevalent, as can be 
observed in most conductors [11]. However, very recently in the history of conductivity 
investigations, superdiffusive and even ballistic motion have been produced in laborato- 
ries. This introduces a new and important field of investigation [12-17]. Indeed, we can 
find reliable reports on ballistic conductivity in carbon nanotubes [12, 13], in semiconduc- 
tors [18], and in semiconductor superlattices with intentionally-correlated disorder [14, 15]. 
We discuss here the violation of MH and EH in ballistic motion. 

For ballistic conduction fl(z) oc cz, where c is a number without dimensions. The limit 
becomes 

lim R(t) = (1 + c)- 1 ^ 0, (9) 
and the MC is violated. This implies that the EH and the FDT are violated [1]. 
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Fig. 1. Plot of the relaxation time r v as a function of the index v for values of cod/ J equal 
to 2, 5, 8, and 10. Note the growth of the maximum as the ratio ujdM increases. 



Briefly, for some processes like anomalous diffusion the Lee condition Eq. for the EH 
is too restrictive and in order to have the EH, we need only the MC. Moreover, the integral 
is finite only for normal diffusion, a = 1. For anomalous diffusion in the range < a < 2, 
the integral is either null (0 < a < 1) or infinite (1 < a < 2). In all these cases, the MC, 
the EH and the FDT are valid. At the extreme limit u = — 1, or a = 0, the system behaves 
as a localized harmonic oscillator and violates the MC, the EH and the FDT; for a = 2, 
ballistic motion, the system also violates MC, and consequently the EH, and the FDT fail. 

Conclusion - In this work we revisited the problem of the validity of the EH and that of the 
MC, and we obtained an agreement between our previous result [1] and recent results using 
the recurrence relations method [4]. However, for anomalous diffusion our condition can 
be less restrictive than that of Lee. Since diffusion is a main phenomenon in physics we use 
that as our starting point. The method used by Lee [2-4] is quite general and applies to any 
response function. In this context, the equivalence between this and our results strengthens 
both. The violation of the EH is exhibited for ballistic and harmonic motions. We have not 
focused deeper on real complex systems; we chose to follow easy-to-understand concepts 
where limits can be analytically obtained. This gave us a good framework for analyzing 
more complex structures. Nonlinear dynamics is a field which deserves much attention; 
in particular, the coalescence of trajectories has been intensively studied in the last few 
years [19-21]. The restriction of the degrees of freedom there may confirm the hierarchy 
exposed here. We also expect that new mathematical methods [22] may bring alternative 
proofs to the problem. 
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